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ABSTRACT 


HYDROFOIL SHIP 
CONSTANT LIFT CONTROL SYSTEM 
PART I: HYDRODYNAMICS 


MICHAEL ROY TERRY 


Submitted to the Departments of 
NAVAL ARCHITECTURE and MARINE ENGINEERING 
and MECHANICAL ENGINEERING on 
23 May 1969, in partial 
fulfillment of the requirements 
for the Degrees of NAVAL ENGINEER, 
and MASTER of SCIENCE in MECHANICAL ENGINEERING, 


This paper analyzes the unsteady hydrodynamics of 
1) FULL INCIDENCE; 2) TRAILING EDGE FLAP; and 3) TAB FOIL 
hydrodynamic control devices, The device is to be used in a 
constant lift, wave alleviation control system on an OPEN- 
OCEAN HYDROFOIL SHIP. Each device 1s compared in POWER, DRAG, 
and WEICHT. The conclusions reached are: 


1. Wave Alleviation is an important aspect of 
hydrofoil ship control system design. 


2. Unsteady Analysis is vital to understandini 
the response of a foil in waves. 

3. Negative Control System Power is required for 
wave alleviation by a tab foi: with a lagging 
tab. 

4, Trailing Edge Flap reauires less overall 
system power than a Full Incidence Foil. 

S. Hydrofoil Ship Operation Policy will be greatly 
ateccel by unsteady hydrodynamics. 

6. Unsteady Foil Response is not a unique function 
of encounter frequency. 


Thesis Supervisor: Damon Ellis Cummings 
Title: Assistant Professor of Naval Architecture 
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SYMBOLS 


NACA camber loading parameter 

a hinge point in % of b 

Foil area 

Ssemichord 

Chord 

B hinge point in % of b 

Designates circulatory 

Lift coefficient L/((1/2)pV7A) 
Lift coefficient at ideal angle of attack 
Theodorsen 2=-D unsteady lift function 
source position 

Dra 

Designates encounter 

Gravitation constant ft*/sec 

Keel to foilborne waterline height 
aT 

Main foil inertia/ft span 

Flap inertia/ft span 

Bessel Function of order n 

Reduced frequency 

Lift/ft span 

Length 

Source strength/ft span 

Moment/ft span 


Mass of the ship 





NC 


Re 


wn BR RE 


Wt We 


Denotes non-circulatory 

Flow rate 

Span 

Time 

Thickness 

Period 

Ship Speed 

Wave Phase Velocity 

Denotes wave 

Distance in horizontal direction 


Longitudinal distance from reference point to 
center of gravity 


Distance in vertical direction 

Distance traveled in vertical direction 
Vertical velocity 

Vertical acceleration 


Vertical force 


Pitch angle of incidence about ba 
Pitch angle of incidence magnitude 
Incidence angular velocity 
Incidence angular acceleration 
Pitch angle of flap about bc 

Pitch angle of flap magnitude 

Flap angular velocity 


Flap angular acceleration 
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Ds 


Ship displacement 
Increment of change 

Small increment 

Propulsive coefficient 
Wave height 

Wave height magnitude 
Wave orbital velocity 

Wave orbital acceleration 
Ship pitch acceleration 
Wave length 

Bis LHMIG. < 

Fluid density in slugs/ft?® 
phase angle 

Lo 

un 

Phase angle between a, Pp 
Angle of attack of ship to waves 
Stream function 


Circular frequency radians/sec 


APL 
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span 


Total hydrodynamic moment about be per unit 
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Plotting symbol for PAVEA 
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Ow 

Plotting, symbol for PAVEB 
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Upwash velocity n 


Wd) 
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Plotting symbol for ,B 


Plotting symbol for ,a 


a hinge line in % of pb 
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I. INTRODUCTION 

Renewed interest in the use of hydrofoil lifting sur- 
faces to raise a ship above the surface of the water 
occurred in the 1950's. Since that time, development of 
complete hydrofoil ship systems has progressed at a rapid 
rate. During the recent development stares, the encounter 
between ships and open-ocean waves has identified the pro- 
blem of ship motion control as being critical. The foil 
encounters lift chanses caused by wave particle orbital 
motion. The wave particle orbital velocity is superposed 
on the uniform in-flow velocity (ship's speed through the 
water). The superposition results in an angle of attack 
change relative to the foil and, therefore, changes its 
lift. The changes in lift are disturbing forces to the 
Ship and occur at the encounter frequency of the ship meeting 


the waves. Figure (1) depicts this situation. 
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figure 1 


This paper will address only open-ocean waves and ocean- 


going hydrofoil ships. 


During; the early stages of development, attention was 





focused on the choice of foil configuration for best results. 
From that lively discussion and experimentation, two major 
design configurations developed. The first is commonly 
referred to as surface piercing, and the second as fully 


Submerged. Figure (2) depicts these configurations. 


surface piercing fully submerged 
netine ws, 2 


iD 





This paper will address only the fully submerged configura- 


tion. Hydrofoil ships with fully submerged foils are inher- 
ently unstable in heave, pitch, and roll when the hull is 
out of the water, because the foil has very little sensiti- 
vity to depth below the water surface. Therefore, a control 
system must be provided. 

While follborne, two modes of operation are possible. 
The first is the platform mode where the ship's motion is 
maintained at constant height above the mean water level, 
at constant attitude about the height, and with no accelera- 
tions. The platform mode is only possible in wave heights 
that do not exceed the keel-to-foilborne waterline distance. 
The second is the contour mode where the ship's motion is 
maintained at constant height above the local water surface 
at two points along the ship's length. The contour mode 
is only necessary in wave heights that exceed the keel-to- 


follborne waterline distance. Figure (3) depicts these 


modes. 
platforming 
— a 
contouring 7 
a) 


LGtise > 
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this) paper will address only the platform mode and the wave 


systems that allow the use of that mode, 
The design of motion control for a hydrofoil ship must 


meet the three basic objectives outlined in Figure (4). 


STABILITY (Maintaining Commanded 
Equilibrium Attitude) 


MANEUVERABILITY (Changing to a New 
Commanded Equilibrium) 


WAVE ALLEVIATION (Reducing the Effect 
on Lift due to Waves) 


figure 4 


Most attention in the past has been concentrated on stability 
and maneuverability control functions. The effect of waves 
on the ship has been treated as a disturbance from equili- 
brium, and the stability elements of the control system have 
been sized to be able to correct the ship's attitude and 
reduce the accelerations to a tolerable level. In contrast, 
this paper will address only the wave alleviation objective. 
The control system is usually composed of elements as 


shown in Figure (5). 
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SENSORS (ELECTRONIC ) 


AUTOPILOT/COMPUTER (ELECTRONIC) 


SERVO (FLUID POWER) 


HYDRODYNAMIC DEVICE (LIFTING SURFACE) 


figure 5 


This paper will address only the Hydrodynamic Device. 
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Ti. PROgsCT 

The objective of the above focus is now described. An 
open-ocean hydrofoil ship will operate, at times, when wave 
systems are present for the complete voyage. Operation in 
waves is, therefore, an important aspect of the design. Would 
the control system perform better than present day systems 


if the approach to its design considered only wave allevia-~ 


tion as the objective? The hypothesis is that the power 
level of the control system operating in waves could be 
Significantly reduced if this approach was taken. Also, the 


wave alleviation control system could be integrated with a 
present day standard stability and maneuverability control 
system to produce a control system with overall performance 
superior to today's design. The stabilizing and maneuvering 
elements would effectively be operating as if they were in 
calm water, when in fact the hydrofoil ship is flying in 


open-ocean waves, Figure 6. 


jr ne 


INflow 
figure 6 


i 





In the process of designing the optimum control system for 
wave alleviation, different solutions must be tested. The 
criteria best suited for this test is minimum overall system 
power. The three elements that contribute directly to over- 


all system power are shown in Figure (7). 


POWER Power to Operate 
the Control Surface 


DRAG Power to Overcome the Drag 
of the Hydrodynamic Device 


Power to Lift, Through the 
WEIGHT Lift/Drag Ratio, the Weight 
of the Control System 


Matinee 7, 


For different designs, any one of the above elements of 
overall control system power may predominate. The task of 
Gesigning the control system can be broken down into the 
design of the four basic elements--sensor, autopilot, servo, 
ana hydrodynamic device as shown in Figure (5). "“Hydrofoil 
Ship Constant-Lift Control System: Part I Hydrodynamics" 


Will address only the fourth element—-liydrodynamic Force 


Producing Device. 
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In summary, Figure (8) shows the steps in focusing this 


thesis within the overall project. 


RIVER INLAND COASTAL 
SEA 

SURFACE 

PIERCING 


CONTOURING 


‘STABILITY MANEUVERABILITY 
SENSORS AUTOP ILOT/ SERVO 
COMPUTER 


figure 3 
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Pol, 


IYDRODYNAMIC DEVICE 
Ae DEVICES 

The first step in the plan to find the minimum 
power device is to develop as many candidates as possi- 
ble for closer scrutiny. In general terms, the hydro- 
dynamic device must effect a force. Although there 
are a number of ways of effecting a force, this paper 
will consider only those that develop lift. This is a 
reasonable specialization because the hydrofoil ship 
utilizes lifting surfaces for lifting itself out of the 
water and for stability and maneuverability. 


Table (1) is an attempt at classifying the candi- 


dates, This paper will pursue only the following 


candidates: 


le. Full Incidence 
2. Trailing Edge Flap 


3. Tab Foil 


Appendix (A) describes the remaining candidates and 
discusses some of their characteristics in qualitative 


terms, 


ee 
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Before proceeding with a quantitative analysis of the 
three candidates, a reference hydrofoil ship should be de- 
fined for use as a common denominator for each design. 
Table (2) lists the parameters of the reference hydrofoil 


Ship chosen, 


TABLE (2) 
Displacement A = 300 Tons 
Length Overall L = 220 Ft 
Design Speed U = 40 kts 
Design Lift/Drag L/D = 7 


Design Propulsive Coefficient 16 = 6 

Each of 2 Main Lifting Foils NACA 16-408 
Thickness/Chord t/co = .06 
Mean Line a= 1.0 


Design Lift Coefficient Cyo = 4 


Chord c = 8,58 ft 
Semichord b = 4,29 ft 
Span s = 21.25 ft 
Wetgeht wt. = 7 tons 


Foil Inertia/Ft Span I, = 140.0 Slug Ft?/ft 
1/4 Chord Trailing Edge Flap 

Chord = 2,145 ft 

Wt. =s 1.75 tons 

Flap Inertia/Ft Span In = 8.80 Slug Ft?/ft 


Keel to Foilborne Waterline Height h = 14.0 ft 





IV. POweR 

The first and most Important component of overall con- 
trol system power is the power to drive the hydrodynamic 
device throuyn motions that will effect platforminy; in 
waves. A reasonable model of motion in the vertical plane 


for the hydrofoil ship is 


Mz + Mx,6 = Z(waves) + Z(Ship + Z(Control CL) 
Motion) Motion) 
the effects of waves, ship motion, and control motion are 
modeled as linearly superimposed. For the platforming mode, 
the equation is reduced to 


Zero = Z(waves) + Z(control (2) 
motion ) 


The objective is to cancel the lift caused by waves by 
generating an egual and opposite lift caused by the control 


surface. From the required motion of the control surface 
and the resultant forces on it, the power required can be 


calculated 


Power = (motion velocity)(summation of forces in (3) 
the direction of the velocity) 
The next task is to look at the effect of waves ona foil. 
A, MQUASI-STEADY LIFT DUE TO WAVES 
From a qualitative, quasi-steady point of view, the 


lift due to a wave passing a fully-submerged hydrofoil is 


eS = 





causeG by the orbital velocity of the water particles super- 
posed with the velocity of the foil through the water. for 
a steady foil, the orbital velocity appears as an angle of 


attack change for the incominz flow, Figure (9). 


lift 


y (“of aa 


MURS = 


For the hydrofoil siip, tie lisiiting wave condition for 
continued platforming ise-total amplitude of the wave equal 
to the keele-to-foilborne waterline height. For the parti- 
cular hydrofoil ship chosen, this means a wave half-height 
of 7.0 feet. For an open-ocean wave, the common wave length 
for this amplitude is 280.0 feet. At the ocean surface, at 
a given point, a reasonable model of the wave vertical posi- 


tion is 
n= Re{n ,exp(-L27V/rA+t )} 


For a deep water wave 


ale 





or 


n= Re{n,exp(-i/2mg7X + t)} (6) 
The wave vertical particle velocity is then 
A = Re{n,(-1/2mg/7i)exp(-i/2mg7X © t)} (7) 


The amplitude of the orbital velocity is then 


In| = no¥2ag7X (3) 
For the chosen wave 
Inj = 5.85 ft/sec (9) 


With the hydrofoil ship's velocity at 67.5 ft/sec, then the 
maximum angle of attack change due to vertical orbital velo- 


eacvy is 
Tan a = 5295/67.5 (10) 


or 


a = 5,.06° (11) 


The maximum inflow velocity change due to horizontal orbital 


velocity is 
U,; = 67.5 + 5.95 ft/sec (12) 


To determine the effect on lift, we must look at the lift 


coefficient Cy. For the 5.06° angle of attack change, flat 
plate model of the foil gives 
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9C,/da a ai (13) 


or 
AC, = 2mda (14) 
AC, = 0546 Crs) 
SO 
C_, = .4 + .546 = .946 (16) 


1 


For the hydrofoil ship chosen, the ratio of the lifts is then 
L,/L, = 2.36 (17) 


In contrast, for the 5.95 ft/sec inflow velocity change, the 
ratio of the lifts is 


Ly) / by (UU (18) 
Ly/Ly = 1.19 (19) 


From the above, the effect on lift due to the vertical compo- 
nent of orbital velocity is far greater than that due to the 


horizontal component of orbital velocity. This paper will 


only consider the effect of the vertical comnonent of orbital 
velocity. 


The use of quasi-steady analysis is not adequate for 
the case of the hydrofoil ship chosen. The encounter 
frequency, semichord, and ship's velocity combine to sive 


a reduced frequency 


k =w,b/U (20) 





in the range where unsteady effects are most important, 
as will be shown in Section (B). Figure (10) describes 


the concept of encounter frequency wW,. 


Wc os 





figure 10 


here 
ena 2n/l, = 27U,/L, 
where 
L, = A/Cosy 
and 
OR (V/Cosp) + U 

SO 

Ong (2n/A)C(V + UCosy) 
but 


V2 = gpr/eon 


therefore 


a f ‘ 
r ome / Ww, 


30 


r 


Ce 


(22) 


(23) 


(24) 


(26) 





or 


2s 2 e S 
We Wi WW UCosw/gZ 


For the wave and ship chosen 


We © 085 + .72(67.5)(1)/32.2 = 2.36 Rad/Sec 


Semichord b = 4,29 ft 
Speed U = 67.5 ft/sec 


sO 


hoe 2. 36(4529)/67 «5 =e 


(27) 


(28) 


a9) 


or k is in the most critical range for unsteady effects as 


shown in Section B. 


Be. UNSTIADY LIFT DUE TO WAVES 


figure 11 


ot 








The following, is an adaptation of Chapter 5 of 
Reference (1) for this particular problem. Please refer 
to Table (3) in this paper for the following discussion, 
A model of a wave passing a steady foil, shown in Figure (11), 
is expressed by Equation (30). The disturbing upwash velo- 
city UPW is equal to a complex phasor with magnitude--|n|, 
phase angle (relative to wave height)--¢,, and being driven 
around the phasor diagram with angular velocity--w,. 

Assuming negligible effects from free surface, negli- 
gible energy radiation on the free surface; idealized 2-D 
flat plate theory, with no separation or cavitation, the 
unsteady lift is expressed as the quasi-steady lift modified 
by the Theodorsen Function and additional Bessel Functions. 
These functions are expressed in Equations (31 to 34). ‘The 


Theodorsen Function is shown in “iLyture (12). 





Im figure 12 


From igure (12), it can be seen that the Theodorsen Function 


reduces in magnitude (from 1 to 1/2) and has its maximum 
lagging phase anple change near k = .153;--exactly where the 


nydrofoil ship is operatins:. 
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The expression for the complex phasor LIFTW is shown 
in Equation (35). The complex phasor moment MOMTW due to 
LIFIW acting about point ba (Figure (11)) is expressed in 
Equation (36). The complex phasor moment. HOMTBW due to lift 
on the flap acting about point be (Pigure (11)) is expressed 
in Equation (37). These functions are best visualized with 


a phasor diagram. Figure (13) shows a phasor diagram scaled 





to illustrate the relationship between UPW, LIFTW, MOMTW, 
and MOMTBW for w. = (+). 
NOTE ; 
1. The diagram rotates clockwise. 
2. The unwash velocity leads the wave height by 90°. 
3. LIFTW lags UPW. 
4, MOMPY is in phase with LIFTW. 


MOMTBW is 180° out of phase with LIFTW. 


Oa 
© 


6. LIFT and MOMTW would have been larger in magnitude 
and in phase with UPW, if k = 0. 

Figure (14) shows a phasor diagram again scaled to illustrate 
the same relationships for w, = (-), 

NOTE: 

1. The diagram rotates counterclockwise. 

2. The upwash velocity lags the wave height by 90°. 

3. Again LIFTW and MOMTW lag UPW. 

4, MOMTBW is 180° out of phase with LIF@W. 


Figure (15) illustrates the meaning of we = (+,-). 
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Cs UNSTEADY LIFT DUE TC aN OSCILLATING JTL 

U b ee = 
Cee ne nN ccm’ “ ba 
oe 





+X + MomrT BO : 
figure 16 ie Ee 


1. Alpha Motion 

The following is again an adaptation of Chapter 5 of 
Reference (1). Please refer to Table (4) in this paper for 
the following discussion. A model of full incidence motion 
(Alpha Motion) (Figure (16)) is expressed in Equations (38 
to 40). It is expressed as a complex phasor with areacce 
angle relative to the wave height--~-¢,, and being driven 
around the phasor diagram with angular velocity~-w,. 

Again, assuming neglivible effects from free surface, 
negligible energy radiation on the free surface; idealized 
2-D flat plate theory, with no separation or cavitation, 
the unsteady complex phasor lift (LIFTA) is expressed in 
non-circulatory velocity and acceleration terms, and circu- 
latory position and velocity terms modified by the Theodorsen 
Function. LirfA is expressed in Equation (41). The unsteady 
complex phasor moment MOMTPA due to LIFTA acting about point 


ba (Figure (16)) 1s expressed in Equation (42). The unsteady 
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complex phasor moment MOMTBA due to lift on the flap acting 


about point be (Figure (16)) is expressed in Equation (43). 


Again, these functions are best visualized with a phasor 


dla ram, 


Figure (17) shows a phasor diagran scaled to 


illustrate the relationship between ALPHA, ALPHADOT, ALPHA- 








DDOT, LIFTA, MOMTA, and MOMTBA for we = (+). 





HOTE: 

lL. The diagram rotates clockwise. 

2. ALPHADOT leads ALPHA by 90°. 

3. ALPHADDOT leads ALPHADOT by 90°. 

4, LIFTA lags ALPHA. 

5. MOMTA and MOMTBA lag LITA. 

6. LIFTA will be in phase with ALPIA for k = 0 


because although LIIFTA has contributions frorn 


ALPHADOT and ALPHADDOT, these are zero for k = 0. 


Figure (18) shows a phasor diagram to illustrate the same 


relationships witn we. = (-). 


HOUT ¢ 


ile 
Qe 
3. 
tI 


Tae diagram rotates counterclockwise. 
ALPHADOT leads ALPHA by 90°. 
ALPHADDOT leads ALPHADOT by 90°. 


LIFTA lags ALPHA. 


2. Beta “Motion 


The following is again an adaptation of Chapter 5 of 


Reference (1). Please refer to Table (5) in this paper for 


1) 
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the following discussion. A model of trailing edge flap 
motion (Beta Motion) (Figure (16)) is expressed in Equations 
(44 to 46). It is expressed in terms of the complex phasor 
ALPHA. The ratio of magnitude between BETA and ALPHA is-~ 
Bo/ag; the phase difference between BETA and ALPHA (+ for 
BETA leading) is--Ad; and the phasor is being driven around 
the phasor diagram with angular velocity—-w,. 

Again, assuming negligible effects from free surface, 
negligible energy radiation on the free surface, idealized 
2~D flat plate theory; with no separation or cavitation, 
the unsteady complex phasor lift (LIFTS) is expressed in 
non-circulatory velocity and acceleration terms, and circula- 
tory position and velocity terms modified by the Theodorsen 
Function. LIFTS is expressed in Equation (47). The unsteady 
complex phasor moment MOMTB due to LIFTB acting about point 
ba (Figure (16)) is expressed in Equation (47). The unsteady 
complex phasor moment MOMTBB due to lift on the flap acting 
about point be (Figure (16)) is expressed in Equation (49). 
Again, these functions are best visualized with a phasor 


diagram. Figure (19) shows a phasor diagram scaled to 





illustrate the relationship between BETA, BETADOT, BETADDOT, 
LIFTS, MOMTE, and MOMTBB for we = (+). 

NOTE: Same as Figure (17). 
Figure (20) shows a phasor diagram to illustrate the same 
relationships with w, = (-). 

NOTE; Same as Figure (18). 
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D. REQUIRED FOIL MOTION 


Mode 1 = ALPHA only (Full Incidence) 


Mode 2 = BETA only (Trailing Edge Flap) 


Mode 3 = ALPHA and BETA (Tab Foil) 


The wave is given. The lift due to the wave can be 


solved explicitly by Section B. However, the required con- 


trol surface motion (ALPHA and/or BETA) is only implicitly 


expressed in Equations (41) and (47) of Section C. Both 


Equations (41) and (47) are complex expressions. Each one 


represents two equations either as the real and imaginary 


parts or the magnitude and phase angle. 


The 


1. 


2 
36 


stratery used to solve for ALPHA and/or BETA is: 
Substitute in Equations (41) and/or (47) expressions 
for ALPITA and/or BETA and their derivatives in terms 
of an unknown ad and $,. (B,)/a, and Ad known con- 
stants ) 

Expand multiple angle expressions. 

Separate the real and imaginary parts of the modi- 
fied Equations (41) and/or (47). (-ReLIFTW, 
~ImLTF iy ) 

Divide out a9, Cosd,,p, and b. (=ReLIFTW/(aoCos¢,pb). 
Group the remaining terms into terms containing 


Tan $1--4, X,, and constant terms--Y,, X,. 





-ReLIPIW/(a,Cos¢,pb) = Tangd,(Y¥i) + (Y2) (50) 
-ImLIFTY/(aoCosdipb) = Tang: (X,) + (X,) (51) 
6. Divide the real equation by the imaginary equation. 


~ReLLFTW/(-ImLiFiW)= R 
= (Tang, (Y, )+(Y,))/(Tand, (X,) + (X,)) (52) 


7. Separate out Tan 4, 
Tan 6, = (Y2—RX2)/(RX,-Y, ) (53) 


8. Solve for 41. 


9. Magnitude a, is 
Gy = -ReLIFTW/[Cosd, pb(Tang,(¥,) + (Y,))] (54) 


10. The problem is solved. 
See Appendix C for details. With ALPHA and BETA, the equa- 


tions in Section C will give lift and moment from the motion. 


E. SUMMATION OF HYDRODYNAMIC FORCES 
Summation of LIFTW, LIFTA and/or LIFls should equal 
zero. Summation of MOMTY, MOMTA and/or MOMTD gives the 





total hydrodynamic moment about ba. Summation of MOMTBW, 
MOMTBA and/or MOMTBB gives the total hydrodynamic moment 


acting on the flap about be. 
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F. REQUIRED APPLIED MOMENTS 

The control system must apply a moment (MOMAP) on the 
control surface to produce the required motion. The dyna- 
mics of this interaction is expressed in Equations (55) 


and (56). 
I,° ALPHADDOT HADDOT = MOMTW + MOMTA + MOMTB + MOMAPA C55) 


I, °BETADDOT = MOMTEW + MOMTDA + MONTH BB + MOMAPB (56) 


From these equations, the required applied moments can he 


found. Figure (21) shows these relationships in phasor form. 
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G. AVERAGE POWER 
The average power (PAVE) required to produce the 


motion is expressed in Equations (57) and (58). 





PAVEA = |ALPHADOT| -|MOMAPA|*POWER FACTOR/2 (57) 
PAVEB =|BETADOT| + |MOMAPB|+* POWER FACTOR/2 (58) 


where the power factor is the cosine of the phase angle 


between the motion velocity and the applied moment. 


H. PROGRAM 

- Two programs were written in Fortran IV for use at 
Meeinformation Processing Center on the "IBM System 360". 
The first program writes out the details of the power 
calculation, and the second writes out a minimum of identi- 
fying information and plots the power results on a "Calcomp 
Plotter", An outline of the program is shown on Figure (22). 

The encounter frequency is calculated using Equation 

(27) as shown in Figure (10). ‘The reduced frequency is 
calculated using Equation (20) but taking the absolute 
magnitude for use with Bessel Functions. The YTheodorsen 
Function is calculated using Equation (31) by calling the 
"IBM Scientific Subroutine Package" routines BESJ and 
BESY. Because of the iterative method of calculation used 
in BESY, the reduced frequency parameter k can not be 
less than 0.05. This is guarded against in the program 
but results in errors in the Theodorsen Function for reduced 


frequency below 0.05. 
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The WAVLIF routine calculates the complex phasor 
forces generated by the wave using the analysis of Section 
B. Tne lift due to waves (TIFIW) is brought into the CONTRL 
routine. Here the complex phasor ALPHA and/or BETA are 
calculated to produce an eaual and opposite lift to cancel 
LIFIW. The procedure used is described in Section D. The 
motion (ALPHA and/or BETA) is brought into the two routines 
LIFA and LIFB to calculate the lift caused by the motion, 1) 
as a check on the results of CONTRL; and 2) in order to 
calculate the moments about ba and bc, using the equations 
of Section C. 

So far, the control surface motion and the hydrodynamic 
moments are known. The hydrodynamic moments due to waves 
and control surface motion are summed as in Section E, Next, 
the required applied moments are calculated as in Section F. 
The average power is calculated using the analysis of Section 
G. Note all results are 2-D per foot of span. An important 
requirement is the compatibility between mode specification 
and the input data values for BA = Bo/ao, and DPHE = Ad. 

For Mode 1 ALPHA motion, 

BA = 0.0 DPHE = 0.0 
For Mode 2 PETA motion, 

BA = 1.0 DPHE = 0.0 
For Mode 3 ALPHA and BETA motion 


DPHE # 0.0 





The plotting routines plot the motion variable magni- 
tude and phase angle, and the average power as functions 
of encounter frequency. One plot per foil, the plotter 
treats each wave run as a point and connects each point 
with a straight line. The choice of wave input data should 
be arranged to give a somewhat continuous variation as was 
done for the series discussed in Section I. 

Appendix D contains a listing of the programs and an 
example of their written output. Examples of the plotted 
output can be found in Section I. 
ieee eS ULTS 

The two pieces of input data required for each run of 
the program are the foil data and the wave data. Three sets 
of calculations were performed for the three candidate foils, 
as described by Table (2): 

Mode 1 (Full Incidence) 

Mode 2 (Trailing Edge Flap) 

Mode 3 (Tab Foil) 

The first set is a series of waves with ship speed con- 
stant, wave length constant, and angle of attack (SI) between 
the ship and the waves varying from "dead ahead" to "dead 
astern", Table (6) displays the wave parameters for each 
run of this set. Figures (23), (24), and (25) are the compu- 
ter plotted results for the three folls run with this series. 
The series can be pictured as a slow 180° turn being executed 
by the ship at constant speed in regular waves of constant 


wave length. 


oa 





RUN # 


Oo ON NU Bw eH 


10 


12 


13 
14 


SHIP ANGLE OF ATTACK TO THE WAVE, VARIATION 
HWen WL= U=U SI=y 
(FT) (a) (FT/SEC) (RADS ) 
7.0 280.0 80.0 0.0 
a0 280.0 80.0 23927 
(6 280.0 80.0 07854 
7.0 280.0 80.0 Teor 
7.0 280.0 80.0 1.3745 
7.0 280.0 80.0 1.5708 
7.0 280.0 80.0 1.7672 
TO 280.0 80.0 1.9635 
Tea 280.0 80.0 2.1599 
7.0 280.0 80.0 2.3562 
7.0 280.0 80.0 2.5526 
Pe 280.0 80.0 2.7489 
7.0 280.0 80.0 2.9453 
7.0 280.0 80.0 3.1416 


TABLE (6) 





PHEW= ow 


(RADS ) 


~1.5708 
~1.5708 
~1.5708 
~1.5708 
-1.5708 
~1.5708 
ol, 5708 
-1.5708 
-1, 5708 
-l, 5708 
-1.5708 
-1.5708 





The second set is a series of waves with ship speed (U) 
varying from 90 ft/sec down to 25 ft/sec, wave length con- 
stant, and angle of attack between the ship and the waves at 
"dead ahead" and "dead astern", Table (7) displays the wave 
parameters for each run of this set. Figures (26), (27), 
and (28) are the computer plotted results for the three foils 
run with this series. The series can be pictured as a slow 
deceleration by the ship in ahead regular waves of constant 
wave length and a slow acceleration by the ship in astern 
regular waves of constant wave length. 

The third set is a series of waves with constant ship 
speed, wave length (WL) varying from 280 ft to 140 ft, and 
angle of attack between the ship and the waves at "dead 
ahead" and "dead astern". Table (8) displays the wave para- 
meters for each run of this set. Figures (29), (30), and 
(31) are the computer plotted results for the three foils 
run with this series, The series can be pictured as a slow 
transit by the ship at constant speed in ahead regular waves 
of reducing wave length and a slow transit in astern regular 
waves of increasing wave length. 

Additional trial and error runs were made in order to 
discover a useful motion ratio By /Q, (BA) and a useful phase 
angle difference Ag (DPHE) for Mode 3 (Tab Foil). The search 
is not complete, but BA = 2.0 and DPHE = ~.3927 radians 


appeared best and were utilized for the Mode 3 runs. 
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SPEED VARIATION 


WL=. 
(FT) 


280.0 
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280.0 
280.0 
280.0 
280.0 
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280.0 
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280.0 
280.0 
280.0 
280.0 
280.0 
280.0 
280.0 
280.0 
280.0 
280.0 
280.0 
280.0 
280.0 


TABLE (7) 


U=U 


(FT/SEC) 


90.0 
85.0 
80.0 
75.0 
70.0 
65.0 
60.0 
55.0 
50.0 
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PHEW= ow 


(RADS ) 


-1.,5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
-1.5708 
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RUN # 


OO CON NU SW DY 


10 
12 


HW=n 
(FT) 
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WL=A 
(FT ) 


140.0 
168.0 
196.0 
224.0 
252.0 
280.0 
280.0 
252.0 
224.0 
196.0 
168.0 


TABLE (8) 


U=U 
(FT/SEC) 


WAVE LENGTH VARIATION 


80.0 
80.0 
80.0 
80.0 
80.0 
80.0 
80.0 
80.0 
80.0 
80.0 
80.0 
80.0 


3.1416 
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figure 29 
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All the results are in error for the range of reduced 
frequency, -—.05<k<+.05 , because of the fix-up procedure 
mentioned in Section H. The error is dramatic because of 
the sensitivity of power to the related phase angles of its 
components. The error is zero at k = 0, and beyond the 
range k =+05. The range for k translates into a differ- 


ent ranze of w, for each of the three wave series. The 


e 
reason is that k is composed of both W and U, which are 
not solely functions of one another, Equation (27). For 


agiven U, w, remains dependent on dA and The following 


e 


are the ranges of Ww, which are in error for: 


Series 1. (SI) Variation 

~ .85<w.<+.85 (Rads/Sec) 
Series 2, (U) Variation 

~.94<w <+.29 (Rads/Sec) 
Series 3. (WL) Variation 


Runs not in the range 


The plotted results are somewhat misleading in a part 
of the W, range, because the plotter plots straight lines 
between points. The following are the ranges in W,. where 


there are significant gaps in plotted points: 


1. (SI) Variation 


None 
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2. (U) Variation 
+.28<w,<+1.41 (Rads /Sec) 
3. (WL) Variation 


~-94<w <t2.64 (Rads/Sec) 


Aside from the fix-up error and taking note of the range for 
the wap in points, the plotted results give a reasonably 
good picture of the power required. 

In the case of Mode 3 ALPHA and BETA Motion (Tab Foil), 
the plotter plots the Beta power as Q= ~ PAVEB. This opera-~ 
tion places the Q result relative to the P result (PAVEA) 
such that the distance between them is the total power re- 
quired by the system. If the distance between Q + P is in 
the positive direction, the total power required is positive. 
if the distance between Q2 + P is in the negative direction, 
the total power required is negative. The above is true 


regardless of where the plots lie on the absolute scale. 
THE FOLLOWING ARE COMMENTS ON EACH FOIL RESULT: 


1. Alpha Motion (Full Incidence) Mode 1 

ae (SI) Variation. Figure (23) shows a decreasing 
average power requirement from "dead ahead" around to 
quartering seas where the encounter frequency goes to zero 
(power goes to zero); and then continuing around to “dead 
astern", the power increases again for we = (-). The maximum 


power occurs at SI = 0 and it is positive. 





Tne magnitude of the Alnha motion is constant, because 
the disturbance is unchanged in magnitude and the ship's 
(lift generating) speed is constant. Note, at wo, = 0 
(SI = 112°), the magnitude of Alpha does not go to zero. 

This is because, depending upon what the phase angle of the 
wave is that the ship is running along (crest, back-face, 
trough, or front-face), Alpha will have to be a constant OQ, 
ma, 0, (-). 

The phase angle of Alpha relative to the wave height 
changes from lagging 90° + e« for aoe (+), to leading 90° - ¢« 


foru, = (-). At wo = 0, the phase angle is 90° from the 


e 
wave height; but, more important, it is 180° out of phase 
with the wave upwash velocity UPW as is expected from quasi- 
steady analysis. The basic characteristic of shifting from 
lagging to leading and being out of phase with UPW at we = 0 


is common to all the foils in all the wave series, 


b. (U) Variation. Figure (26) shows increasing 
magnitude for the required Alpha motion, as speed drops from 
90 ft/sec in ahead seas (lesser wig (+)). The magnitude 
decreases as the ship speed increases in astern seas (greater 
We = (-)). 

Tne phase angle goes through a similar variation as in 
the SI series, 

The average power decreases as the ship slows in ahead 
seas. In astern seas, the power dips down through zero at 


We = O and then climbs back up as speed increases. The gap 


76 





between points at + 25 ft/sec is necessary because of the 
limited capability of the lifting foils to support the 

ship at low speed. This is evident by the rapidly increasing 
motion magnitude (shown in Figure (26)) required to produce 


the lift at low speeds. 


c. (WL) Variation. Figure (29) shows the magnitude 
of Alpha motion increases slightly as the wave length 
increases in ahead seas. This is because the increased 
wave emeeh presents an increasing disturbance, The same 
variation appears in astern seas. 

The phase angle changes in a similar fashion as in the 
Bb. series, 

The average power decreases as wave length increases 
in ahead seas (lesser uw, = (+)). In astern seas, the power 
again decreases as wave length increases (lesser we. = (-)). 
Note, there is no effect of C(k) error in this plot, because 


Wo is out of the range, 


2. Beta Motion (Trailing Edge Flap) Mode 2 

a. (SI) Variation. Figure (24) shows the magnitude 
of the Beta motion constant as in the Alpha case. The magni- 
tude of Beta is larger than the Alpha case which is to be 
expected for a 1/4 chord flap compared to full incidence. 

the phase angle goes through a similar variation as in 


the Alpha case, 


ek 





b. (U) Variation. Figure (27) shows the magnitude 
increasing in a similar manner as in the Alpha case. Again, 
the magnitude of Beta is larger than that for Alpha. 

The phase angle change is similar to the (SI) variation. 

The average power follows a similar trend as in the 
Alpha case; but, again, it's sensitivity to error in C{k), 
for the range ~.94<w <t+.29, is much greater. The only point 
that is correct in this range is at w, = 0. Again, it is 
speculated that, if correct values for C(k) were used, the 
power curve would fair into zero at w., = 0 and climb back 


© 
up again in a similar fashion to the Alpha case, 


ce. (WL) Variation. Figure (30) shows the magnitude 
of Beta increasing ina similar fashion as in the Alpha case, 

The phase angle changes in the same manner as in the 
Alpha case, 

The average power increases for increasing wave length 
in ahead seas (lesser we = (+)). This result is in contrast 
with the Alpha result where power decreased for increasing 
wave length in ahead seas. The average power in astern seas 
decreases for increasing wave length (lesser we = (=)). 

Note, there is no effect of C(k) error in this plot, because 


We, is out of the range. 


3. Alpha and Beta Motion (Tab Foil) Mode 3 
a. (SI) Variation. Figure (25) shows the magnitude 


of Alpha as a constant straight line as in the Alpha only 
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case. The associated Beta magnitude is two times the Alpha 
magnitude; however, it is not plotted. This is true for the 
(U) variation and (WL) variation as well. 

The Alpha phase angle has a similar trend, as before, 
but the w, # (+) variations are displaced away from 90°. 
The Alpha motion now lags the wave height by 90° - €« for 
w, = (+), and leads by 90° + € for w= (~). The ramp-type 


e 


Jog near w, = 0 is caused by the plotter moving from w, = 


e 
+.16 + ~.14, The jog would approach a step type for points 
approaching wo, = 0. The Beta phase angle lags the Alpha 


phase angle by 22 1/2° for both w. = (1); however, it is not 


e 
plotted. This is true for the (U) variation and (WL) varia- 
tion as well. 

The results for power are important. As was pointed out 
earlier in this section; if the distance hetween Q (Beta 
power) + P (Alpha power) is in the negative direction, the 
total power to drive the motion is negative (the control 
system must absorb power from the motion). Figure (25)'s 
plot of 9, P is exactly this case. Figure (32) demonstrates 
in qualitative form the relative position of the motion 


phasors and their lift phasors. 
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Figures (33) and (34) demonstrate in qualitative form how 


these phasors contribute to negative power. 
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b. (U) Variation. Figure (28) shows the magnitude 
of Alpha increasing as in the Alpha only case. 

The Alpha phase angle variation is similar to the 
Mode 3 SI case, 

The average total power is negative over the range of 
uw 


except near w. = +1.50. Here the power required by Alpha 


e e 


is positive, but the power required by Beta is just enough 


negative for the total power to be zero. 


c. (WL) Variation. Figure (31) shows a slightly 
increasing Alpha magnitude as in the Alpha only case. 
The Alpha phase angle is changing in a manner similar to 


the SI case, 
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The total average power starts out positive for small 


wave lengths in ahead seas. At w, = +4,28, the total power 


e 
is zero and then goes negative for decreasing w,. The 
average negative power for Beta is not enough negative to 
cancel the positive Alpha power for w, +4.28, Note, the 
C€(k) error does not effect this plot. 

From the above results, it is clear that an under- 
standing of the unsteady effects on foils is not only 
necessary but also rewarding. Table (9) lists some of the 


important power results. 


TABLE (9) 
(Full Incidence) (Trailing Fdge Flap) (Tab Foil) 

Alpha Motion Beta Motion Alpha and Beta Motion 
Maximum + Maximum + Maximum t 
Average Power Average Power Average Power 
ft-lb/sec ft-lb/sec ft-lb/sec 
___per ft span per ft span per ft span 

U + 774.5 + 534.7 ~ 153.5 
WL +1905. 2 +491.2 + 23.6, - 133.8 


S2 





V. DRAG 
The components of drag for a foil are: 

Skin Friction 

Pressure Drag 

Induced Drag 

Wave Drag 

Interference Drag 
Each of these has a contribution from: 

Steady Motion 

Waves 

Foil Motion 
The task of sorting out all of these parts is formidable. 
For a 2=-)D analysis, the induced parts can be eliminated. 
For a small strut, the interference drag can be neglected. 
For a well-submerged foil, the wave drag can be neglected, 
Assuming sinall motions with no separation, cavitation or 
ventilation, the pressure drag can be neglected. 

This leaves the effects of steady motion, waves, and 
foil motion on skin-friction drag. The temptation is to 
linearly superpose these three effects as was done for lift. 
Yet the model for skin friction is the boundary layer. 

Reference (2) describes the effect on foileskin friction 
by a uniform inflow with a harmonic perturbation superim- 
posed. This is a model of the effect of waves on the skin 
friction of a steady foil. The shearing stress leads the 
velocity perturbation (horizontal orbital velocity for the 


wave), and the phase angle approaches 45° for k > », 
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An oscillating cylinder in a fluid at rest is discussed, 
and an oscillating infinite flat plate ina fluid at rest is 
discussed. The problem of a finite chord oscillating flat- 
plate foil in uniform flow is not discussed. No mention is 
made of the advisability of superposing such a solution on 
the effects due to waves. 

Although the required motion of the control surface 
results in zero total lift, the lift distribution is differ- 
ent for the wave and for the motion. This is evidenced by 
the different values for wave moments and motion moments. 
Therefore, the pressure distribution will oscillate back 
and forth across the foil. It is speculated that the problem 
can no longer be solved by linear superposition, even if 
there was a solution to the foil motion effects. The solu- 
tion of the coupled effects of waves and motion appears 
beyond analytical solution today. 

All that is left is skin-friction drag from steady 
motion in uniform flow. The result for the three candidates, 
with the same overall foil geometry, is that there is no 
Significant difference in sidiemerdecton drag. 

This statement is adequate for comparing the overall 
power of the three candidates, However, if comparisons of 
these foils were being made with other systems, such as 
spoilers, the above rationalization is not adequate because 
ventilation and separation play a much different part in the 


drag of such systems. 
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Va. WEIGHT 

Any statement about weight without a first iteration 
of the complete control system design must be qualitative. 
However, for the effect of weight to influence the choice 
of control system, an absolute figure must be found. For 
the purpose of emphasizing that power must be supplied to 
lift the weight of the control system, an approximate 


analysis will be pursued, 


A. REFERENCE FULL-INCIDENCE SYSTEM 
To use as a reference system, the weight of a full-~ 


incidence (Alpha motion) control system is estimated as 


follows: 
Group I . Foil Weight =15,700 lbs = 740 lbs/ft 
Group II Mechanical Drive Links = 650 lbs = 30 lbs/ft 


Hydraulic Servo Cylinders = 650 lbs = 30 lbs/ft 
Hydraulic Fluids and Pipes = 1,070 lbs = 50 ilbs/ft 
400 lbs = 20 lbs/ft 


Hydraulic Pumps and Drive 


200 lbs = 10 lbs/ft 


Group III Autopilot and Computer 


Group IV Sensors a 200 lbs = 10 lbs/ft 


B. COMPARATIVE WEIGHT 
The Group I, III, and IV weights for the three candidates 


will not be significantly different. The Group II weights 
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will approximately vary linearly with average power. For 


comparison, from Table (9): 


Trailing Edge Flap Max Power. ., _bg4.7 . 4g 
Full Incidence Max Power HIGS se 
Tab Foil Max Power | = -154. ey 
Full incidence Max Power 05.2 : 
(Here the maximum negative power is taken to size the 
servo machinery even though it must absorb the power.) 
C. TRAILING EDGE FLAP SYSTEM WEIGHT 
Group II = 70 lbs/ft 


D. TAB FOIL SYSTEM WEIGHT 
Group II = 20 lbs/ft 


Ee POWER 
The weight becomes power through the .lift/drag ratio 
and the speed of the hydrofoil ship: 


Power = (Drag/Lift)*Weight Speed (59) 


For comparative purposes, the lift/drag ratio is 7 and the 
speed is 80 ft/sec. Power results due to weight are as 


follows: 


Full Incidence * 1500 fte-lb/sec per ft span 
Trailing Edge Flap -« 800 ft-lb/sec per ft span 


Tab Foil et 230 ft-lb/sec per ft span 
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VII. CONCLUSIONS AND RECOMMENDATIONS 

In summary, the contributions of power, power due to 
drag, and power due to weight have been developed for three 
candidate hydrodynamic devices for use in a constant lift, 


wave alleviation, control system for a hydrofoil ship. 


A. OVERALL CONTROL SYSTEM POWER 

A servo system efficiency is required to convert the 
control surface power into input power to the system. An 
efficiency of .6 1s assumed and has been applied to the 
maximum power results in Table (9). 


Table (10) shows the results for overall control system 





power. 
TABLE (10) 
Trailin 
Full Incidence Eage Fla Tab Foil 
-ib/sec ftelb/sec ft-lb/sec 
per ft span per ft span per ft span 

Power + 1840, + 890. + 4o, 
Drag 
Power 
Weight 
Power + 1500. + 800. + 230. 
Total + 3340. + 1690. + 270. 
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From this analysis, the tab foil looks very promising for 
use as the hydrodynamic device in a constant lift control 
system. Development of the sensor, servo, and autopilot 
system designs for use with the tab foil is recommended. 
If the choice is restricted to either full incidence 
or trailing edge flap, the trailing edge flap looks most 
promising. The development of a control system utilizing 
a full incidence foil for stability and maneuverability, 


and a trailing edge flap for wave alleviation, is recommend- 


ed. 
But the most important recommendations are: 
1. Wave alleviation should be an important considera- 
tion in the design of an open=-ocean hydrofoil ship. 
2. Unsteady lift analysis should be used during the 
design of the control system for an open-ocean 
hydrofoil ship. 
B. SENSORS 


The control system must be able to sense the magnitude 
of the upwash velocity, the direction, and the encounter 
frequency of the wave the foil is about to see, enough in 
advance, in order to control proper foil motion. Wave height 
information is not enough. It 18 suggested that wave height 


and orbital velocity be measured ahead and behind the ship. 


C. MOTION 
Figure (35) shows full incidence, trailing edge flap, 


and tab foil motion. The amplitudes have been exaggerated, 
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but their relative magnitudes have been maintained. The 
most interesting motion is that of the tab foil. With 

22 1/2° lagging phase angle for Beta (-A¢), the Alpha motion 
is caused by the moments generated by the Beta motion. The 
motion appears "natural" in contrast to the 22 1/2° leading 
system shown below it. Figure (36) shows motions for 90° 
lagging, 180° out of phase, and 90° leading. 


——— (~ wet ) 
~ - 90° oO a 20°AP 


“360° ~270 














\ a ae 
figure 36 
Of these, only the 90° lagging appears "natural", Preliminary 


checks with the computer program for 80 ft/sec speed in 


limiting ahead seas, indicate: 


90° Lagging } Alpha Power << 0 


Beta Power << 0 
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180° Out of Phase } QO << Alpha Power 


Beta Power < 0 


90° Leading } 0 << Alpha Power 


0 << Beta Power 


The "natural" looking, motions require negative power. The 
other motions require positive power. References (3) and 
(4) discuss the use of a tab to initiate aileron motion. The 
tab is moved in the opposite direction from the direction the 
aileron will move. At first, this sounds like the aileron 
tab should be 180° out of phase. However, the Beta motion 
must lead (in time) to cause the motion. This lead (in time) 
in the opposite direction is a lag relative to the aileron 
motion. 

For the hydrofoil ship, lags near 90° result in large 
negative powers to be absorbed (like the flutter situation). 
If small lags around 22 1/2° are used, small negative powers 


result. 


D. EFFECTS ON HYDROFOIL SHIP OPERATION 

In the limiting sea (ahead or astern), slowing down the 
Ship will reduce the control system power but will increase 
the magnitude of the control surface motion. Since most 
foil arrangements have a mechanical limit to foil travel, 
Slowing down may not be the rest operational change. Running 
with the wave requires no foil motion and no control system 


power. Running with the crest or trough requires no control 
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surface angle. Running with the face or back requires full 
- or + control surface angle. Running with the crest is 
recommended because no control surface angle is required and 
maximum submergence of the foils is obtained, 

The angle of attack to the sea will reduce the control 
system power but not the foil motion magnitude. When running 
along the wave, the face or back of the wave present no span- 
wise flow; but the control surface will be at full = or + 
angle. Running along the crest or trough will present a 
small spanwise flow to the foils and struts, but no control 
surface angle will be required. Running along the trough is 
recommended because no control surface angle is required and 
maximum submergence of the foils, on the beam, is obtained. 

Finding smaller waves decreases control surface motion 
in all cases, and will reduce the control system power for 
a trailing edge flap system in ahead seas. But, smaller 
waves increase control system power for the full incidence 
system in both ahead and astern seas. Smaller waves increase 
control system power for the trailing edge flap system in 
astern seas, 

It is recommended that wave alleviation unsteady effects 


be analyzed before determining operational policy. 


EB. STATISTICAL ANALYSIS 
The foil responds to wave upwash velocity as an input. 
Therefore, an orbital velocity spectrum with direction infor- 


mation is required. This is a 3=dimensional-surface from 
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statistical wave data. Correspondingly, the foil response 
spectrum is no longer a single curve but is a 3=dimensional-= 
surface for a given ship's speed. The complete picture is a 
series of response surfaces for a range of speeds, 

It is recommended that an attempt be made at bringing 
the unsteady foil response analysis into a statistical for- 
mat. From this, a complete and realistic description of 


hydrofoil ship energy requirements could be made. 
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APPENDIX A: HYDRODYNAMIC DEVICES FOR SECTION III-A 


I. ANGLE OF ATTACK CHANGERS 
A. TRANSLATING FOIL (SUSPENDED FOIL) 

The foil vertical velocity must match that of the incom- 
ing wave. If the speed of the foil is zero, W, * Wye The 
vertical velocity of the incoming wave will be |n| = nw,. 

The velocity of the foil must be |Z] = |n|; but, if Oo 
than {z| = [n|. In other words, the foil motion would exact- 
ly follow the wave surface. If speed is increased, the 


magnitude of w, increases and |z| < |n|. Such large excur- 


e 
sions would probably not be feasible. A patent for a similar 


scheme is in Reference (5). 


II. CAMBER CHANGERS 
A. PARTIAL CHORD FLUID (JET FLAP) 

Much work has been done on jet flaps for high lift on 
aircraft. The feasibility of pumping water through the foil 
has not been determined as yet. Fluidic controls could be 
employed to change the direction of the jet forvwave allevia- 


tion. See References (6) and (7). 


B. UPPER SURFACE MECHANICAL (MEMBRANE ) 

This scheme would employ a second skin for the upper 
surface of the foil. The space between could be pressurized 
to change the membrane shape. The upper surface lift load 
would have to be taken in tension in the membrane. The 


structure appears to be the greatest problem, 
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C. UPPER SURFACE FLUID (FLUID CAMBER) 

The concept of using a layer of fluid on the upper 
surface to change the apparent thickness to the oncoming 
flow was pursued to some detail early in the project. The 
following is an approximate analysis for such a scheme. 

In order to change the lift coefficient of the foil 
from .4 to .9, the camber must change from 2.2% of chord to 
5.2% of chord, or a change of 3% of chord. In order to 
change the camber by 3%, the thickness must be changed by 6% 
of chord. 

The proposed scheme of using "fluid camber" to change a 
foil's lift characteristic calls for pumping fluid from the 
leading edge on one surface of the foil and retrieving the 


same fluid back into the foil at the trailing edge. 





After some fruitless searching for a potential flow 
model for the concentrated jet at the leading edge, a half- 
plane source and sink were chosen as a simple model of both 
entry and exit of the flow. The stream function for a source 


and sink of equal strength m and separated by a distance d is: 


, “ -2y(1 ) 
y = so Tan ‘Ss ty yt ss (1) 
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The stream function for a uniform flow in the x direction 


at speed Ue is: 
y = Uey (2) 


For the superposition of these flows: 


z m ~1 -2y(1/2)d 
vr vey + Beten” sresectzeyaye (3) 


We need the half plane thickness of the ellipse formed at 
» = 0 to be 6% of chord, or y/d = .06. With p= 0, x = 0, 


y = .06d, we can solve for Mm. 
m = a ~06d)2T 
Les Tanwi(Tan B) (4) 
- ~2(.06)(1/2)d 
zen) Cees (5) 


where 


Tan(B) = .243 (6) 
B = -166.67° rH (7) 
m= ,1295+deU, (8) 


But we need only 1/2 this flow but across the span of the 


foil. 
B= = .0648* Useds (9) 
Qs .0648°U,¢A ft®/sec (10) 
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For an ocean hydrofoil ship of 300 tons at 40 knots 
with lift coefficient of .4, requires this area: 


ae L | (11) 
(T72Ypv2C, 


5 00° 2240 s 2 
: COSTS CCT BYICT) BEORRE (12) 


This gives us a required flow rate for the "fluid camber" 


system of: 


9 = .0648 (67.5)(368) = 1610 ft*/sec 
or = 7.22 x 10° gal/min (13) 


The "fluid camber" scheme is not feasible for the full scale 


hydrofoil ship. 


III. CAMBER AND ANGLE OF ATTACK CHANGERS 
A. ROTATING FOIL/PARTIAL CHORD FLUID (JET TAB) 

This scheme would utilize the jet momentum and effect on 
the overall flow to force the Alpha motion. Reference (7) 
discusses such a scheme. Hopefully, the flow rate required 


in the jet could be reduced from that ina jet flap system. 


Be TRANSLATING FOIL/PARTIAL CHORD MECHANICAL (SUSPENDED FLAP) 

C. TRANSLATING FOIL/PARTIAL CHORD FLUID (SUSPENDED JET FLAP) 
Both schemes would use relative motion of a part of the 

chord to reduce the excursion magnitude of the foil. These 


schemes seem worth further analysis. 
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IV. LIFT SPOILERS 


A. MECHANICAL (SPOILERS ) 

Mechanical Spoilers have beén used in aircraft as 
ailerons. The scheme for a hydrofoil ship wave alleviation 
would require the foil to be at the angle of attack necessary 
to produce the maximum lift needed. Spoilers would then 
reduce that lift down to the minimum lift needed. The power 
and weight for such a system will be small, but the drag is 


expected to be great. 


Be FLUID (VENTILATION) 

The scheme of ventilating portions of the lift distribu- 
tion has been developed by SUPRAMAR. However, the same pro- 
blems of drag enter as in Section A. It may be that overall 
system power for such a scheme would be small compared to 


full incidence, 
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APPENDIX B 


DERIVATION OF MOMTBW 
FOR SECTION IV~B 
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APPENDIX C 


DETAILS OF THE IMPLICIT SOLUTION 
FOR ALFA AND BETA FOR SECTION IV-D 
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APPENDIX D 


PROGRAM LISTINGS FOR SECTION IV=H 
I, COMPUTER PROGRAM FOILDYN (WRITING) 
II, COMPUTER PROGRAM FOILDYN (PLOTTING) 
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